Algebraic Structure of Lepton and Quark Flavor Invariants and CP Violation 
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Lepton and quark flavor invariants are studied, both in the Standard Model with a dimension five 
Majorana neutrino mass operator, and in the seesaw model. The ring of invariants in the lepton 
sector is highly non-trivial, with non-linear relations among the basic invariants. The invariants 
are classified for the Standard Model with two and three generations, and for the seesaw model 
with two generations, and the Hilbert series is computed. The seesaw model with three generations 
proved computationally too difficult for a complete solution. We give an invariant definition of the 
CP- violating angle & in the electroweak sector. 
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I. INTRODUCTION 

The observation of neutrino oscillations requires that 
the Standard Model [1-3] be modified to account for neu- 
trino masses. The leading theory of neutrino mass is 
the seesaw model [4] , which contains additional fermions 
which are singlets under the SU(3) x SU (2) X U(l) gauge 
group. An attractive feature of the seesaw theory is 
that it naturally gives rise to leptogenesis [5] in CP- 
violating heavy neutrino decay. The generated lepton 
asymmetry then produces a baryon asymmetry via Stan- 
dard Model sphaleron processes. Interestingly, the light 
neutrino masses favored by experiment are precisely in 
the range needed to produce a baryon asymmetry of the 
right magnitude [6, 7]. The baryon asymmetry is known 
to 10% accuracy from the precision cosmic microwave 
background data of WMAP [8]. 

At energies well below the mass scale M of the heavy 
fermion singlets in the seesaw model, one constructs a 
low-energy effective theory obtained by integrating out 
the heavy Majorana neutrino singlets. The Lagrangian of 
the low-energy effective theory is given by the renormal- 
izablc Lagrangian of the Standard Model plus additional 
higher-dimensional terms obtained from integrating out 
the heavy neutrinos. The leading term is a dimension- 
five operator [9] which produces a Majorana mass term 
for the neutrinos of the Standard Model when the Higgs 
doublet acquires a vacuum expectation value v. This op- 
erator is the unique dimension-five operator which can 
be constructed from Standard Model fields. Thus, it is 
natural for this dimension-five operator to be the first ob- 
served effect of new physics beyond the Standard Model. 
The low-energy effective theory contains additional op- 
erators at dimension six [10, 11]. The leading effect 
of these operators is a flavor-nondiagonal correction to 
the weakly-interacting neutrino kinetic energy term after 
electroweak symmetry breakdown. This contribution re- 
sults in a small 0(v/M) nonunitary contribution to the 
lepton mixing matrix C/pmns- Unfortunately, for GUT- 
scale values of the seesaw scale M, this nonunitarity of 
f^PMNS is far too small to be observed experimentally. 

For the purposes of this paper, the Standard Model 
low-energy effective theory is the SU(3) x SU{2) x U(l) 



gauge theory with only left-handed doublet neutrinos, 
plus an additional dimension-five gauge invariant oper- 
ator which gives a Majorana mass to the neutrinos af- 
ter spontaneous symmetry breaking, and the high-energy 
(renormalizablc) theory is the seesaw model. 

Flavor violation of quarks and leptons by Standard 
Model weak interactions is parameterized by unitary 3x3 
matrices, the CKM matrix in the quark sector and the 
PMNS matrix in the lepton sector. The fundamental pa- 
rameters in the Standard Model are the quark and lep- 
ton Yukawa coupling matrices and the flavor matrix of 
the dimension- five Majorana mass operator. The fermion 
masses and mixing angles are derived quantities, ob- 
tained from the eigenvalues and eigenvectors of the flavor 
matrices in the low-energy theory. In the mass eigenstate 
basis, one still has the freedom to make phase rotations 
on the fermions fields, which leads to the redefinition of 
the CKM matrix 



V 
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in the quark sector, where <&[/ = diag(</> u , </> c , </> t ) and 
$£> = diag(0d, 4>si 4>b)- Physical quantities are basis 
independent, and must be invariant under the rephas- 
ing Eq. (1). CKM rephasing invariants have been 
studied extensively in the literature [12-15], the best- 
known example being the CP-odd Jarlskog invariant 
J = Im V\\V2iy iny^i- Rephasing invariance also exists 
for the lepton mixing matrix. In a previous paper [15], 
we extended the analysis of rephasing invariants to give 
a complete classification of these invariants for the Stan- 
dard Model, and for the seesaw model. 

The parameterization of the flavor structure in terms 
of masses and mixing angles is convenient for computing 
decay rates and scattering amplitudes. However, if one 
wants to understand the origin of flavor structure, the 
more fundamental quantities are the flavor matrices in 
the Lagrangian from which the masses and mixing angles 
are derived by diagonalization. A well-known difficulty 
is that the flavor matrices are basis-dependent, since one 
can make unitary transformations on the quark and lep- 
ton fields in the Lagrangian. For example, the Yukawa 
matrix for charge 2/3 quarks transforms as 



Yu 



Ur. 



Yu U Q 



(2) 
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where Uq and U^c are unitary transformations on the 
quark doublet and singlet fields. One cannot directly 
compare a mass- matrix prediction with experiment, since 
the mass matrices are basis-dependent. Observable quan- 
tities must be independent of this change of basis, i.e. in- 
variant under Eq. (2), and such quantities are sometimes 
referred to as weak basis invariants [16, 17]. One can 
check the predictions of a flavor model by comparing in- 
variant quantities with their corresponding experimental 
values. 

Classifying invariants also is important in analyzing 
theories which explain flavor by a dynamical mechanism. 
The idea can be illustrated by a simple example — con- 
sider a low-energy theory which has a 3 x 3 hermitian 
traceless flavor matrix X which transforms as an SU (3) 
adjoint, X -> UXUK Ima gine that X is a dynamical 
variable in some high-energy theory, and that the low- 
energy value of X is given by minimizing an effective po- 
tential V(X) generated by the high-energy theory, ft is 
well-known (see Sec. V) that the only independent invari- 
ants are Ii = (X 2 ) and ^3 = (X 3 ), where (*) denotes a 
matrix trace, so the potential can be written as V(l2, 13) 
and minimizing it leads to the equation 



= 



dV 



(3) 



where X = X a T a . Eq. (3) implies that X a = 
k d a bcX b X c where the constant of proportionality is k = 
-(3/4)(dV/dI 3 )/(dV/dI 2 ) evaluated at the minimum. 
The solutions of this equation are either (i) the trivial 
solution X = 0, or (ii) X can be brought to the diagonal 
form 



X 




V3 
k 



(4) 



with symmetry breaking in the T 8 direction. Thus the 
SU(3) symmetry is cither (i) unbroken or (ii) broken to 
the SU (2) x U(l) subgroup. Symmetry breaking to U (l) 3 
is not allowed. Examples of this type were studied in the 
early literature on unified theories [18, 19] in the context 
of understanding gauge symmetry breaking patterns by 
minimizing Higgs potentials. A recent example from fla- 
vor physics needing the classification of invariants can be 
found in Ref. [20]. 

There is an extensive literature on quark and lepton in- 
variants (see, e.g. [16, 17, 21-24]). The main emphasis in 
previous work has been the study of CP violation. CP- 
violating invariants analogous to the Jarlskog invariant 
were written down. The vanishing of the CP- violating 
invariants was sufficient to guarantee the vanishing of CP 
violation in the CKM and PMNS mixing matrices. 

In this paper, we take a different approach, studying 
all the invariants, and treating the problem using the 
methods of invariant theory [25-27], which considers the 
ring of polynomials that are invariant under the action 
of a group. Polynomial invariants also are the relevant 



objects for physics applications, since an effective La- 
grangian is written as a polynomial in the basic variables 
which describe the theory. 1 A basic result of invariant 
theory is that the ring of invariants has a finite number 
of generators. There can be non-trivial relations among 
the invariants, known as syzygies [28], so that the invari- 
ant ring need not be a free ring. The number of invariants 
of a given degree is encoded in the Hilbert series. The 
complete classification of the invariant ring is, in general, 
a very difficult computational problem. 

In this paper, we study the invariants of the Standard 
Model low-energy theory and the seesaw theory in both 
the quark and lepton sectors. In the quark sector, the 
complete structure of the invariant ring is given, and the 
relation between the polynomial invariants and rephasing 
invariants also is given. The structure of the invariant 
ring in the lepton sector is considerably more involved 
than in the quark sector. The classification of lepton 
invariants is given for the low-energy effective Standard 
Model theory for two and three generations. For the 
high-energy seesaw theory, the classification is given for 
two generations. For three generations, we have been 
unable to completely classify all the relations or to de- 
termine the Hilbert series because the problem is com- 
putationally too difficult. The simpler invariants (i.e. of 
small degree) are given for this case. 

The paper is organized as follows. Section II defines 
the high-energy seesaw theory and its low-energy effec- 
tive theory. The flavor-symmetry breaking matrices and 
i9-angles of each theory are given, together with their 
transformation properties under flavor symmetry and 
CP. Section III defines the mass and mixing matrices 
of the high-energy and low-energy theories. The high- 
energy theory contains three mixing matrices, the quark 
CKM mixing matrix Vckm, its analogous lepton mixing 
matrix V and a mixing matrix for the heavy neutrino 
singlets W . The low-energy theory contains two mixing 
matrices, the quark CKM mixing matrix Vckm and the 
lepton PMNS matrix C/pmns- Section III explains the 
counting of mixing angles and phases for the mixing ma- 
trices for arbitrary numbers of Standard Model fcrmions 
and neutrino singlets in both the high- and low-energy 
theories. Finally, rephasing invariance of the mixing ma- 
trices is discussed. Section IV provides a brief introduc- 
tion to the mathematics of invariant theory that we need 
for our analysis. Several model theories are considered 
to elucidate the mathematical results. The next two sec- 
tions consider the classification of flavor invariants for 
the high-energy seesaw theory and its low-energy effec- 
tive theory. Section V reviews the classification of the 
quark mass matrix invariants, which are identical to the 
quark invariants of the Standard Model for both theo- 
ries. Sections VI and VII consider the classification of 



1 For example, the chiral Lagrangian is a polynomial in the quark 
mass matrix M. 
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lcpton mass matrix invariants for two and three genera- 
tions of fermions, respectively, in both the low-energy ef- 
fective theory and the seesaw theory. The complete clas- 
sification is given for the low-energy effective theory for 
two and three generations. The lepton invariant analysis 
of the full seesaw theory is significantly more complex. 
The complete classification is given for two generations 
of fermions, and partial results for three generations are 
given. 



II. FLAVOR SYMMETRIES 

We consider the SU{3) x SU(2) x U{1) gauge the- 
ory with rig generations of Standard Model fermions 
and n' g generations of gauge singlet fermions (neutrino 
singlets). The fermion multiplets are Qi = (3, 2) 1 / e , 
U? = (3,l)_ 2/3 , D1 = (3,l)i/3, Li = (l,2)_ 1/2 and 
E? = (l,l)i, i = l,...,n g , and Nf = (1,1) , / = 
All fermion multiplets are left-handed Weyl 

n g have a natu- 



1,..., _ 

fields. The fermion multiplets with n' g 
ral embedding in the 16 spinor representation of SO(10), 
so the usual choice is n' g — n g . Theories with n' g 7^ n g 
also are possible, however. Experimentally, we know that 
n g = 3, but there is no experimental limit on n' g . Big- 
bang nucleosynthesis constrains the number of neutrino 
flavors to be less than four; however, this only constrains 
neutrinos which are light enough to be present at tem- 
peratures of order an MeV. 

The flavor symmetry of the fermion sector of the high- 
energy theory is SU (n g ) 5 x U[n' ) x U(l) 2 , since there is a 
separate SU(n g ) flavor symmetry for each of the five mul- 
tiplets Q, U c , D c , L and E c , a U(n' g ) flavor symmetry for 
the singlets N c , and two additional non-anomalous U(l) 
flavor symmetries. Out of the six possible U(l) symme- 
tries, only three linear combinations are non-anomalous 
under SU(3) x SU(2) x U(l): N c number which is in- 
cluded in U(n' g ), (B — L), and (E c + D c - U c ) num- 
ber. The three additional anomalous U(l) groups can be 
treated as symmetries if the three ^-angles 2 $3,2,1 of the 
£[7(3), SU{2) and U{1) gauge groups transform under 
arbitrary chiral phase transformations tp — * e lai "ip on the 
fields ib = Q, U c , D c , L and E c as 



$3 
2 



§3 - n g (2a Q + ajjc 4 
•&2 - n g (3a Q + a L ) , 
1 4 



OLD' 



(5) 



- 6 * Q + ~ 3 au. 



-a L + a E <- 



Eq. (5) does not depend on n' g or ajyc, since N c are gauge 
singlets. With the transformation Eq. (5), the chiral fla- 



2 The i9 angles multiplying FF terms are not to be confused with 
angles 8 of the quark and lepton mixing matrices. There are no 
instantons in the U(l) sector, but the 1? angle can have physical 
consquences in the presence of topological defects. 



vor symmetry becomes U(n g ) 5 x U(n' g ), with a separate 
flavor factor for each of the six fermion multiplets. 

The U(n g ) 5 x U(n' g ) flavor symmetry of the fermion 
and gauge kinetic energy terms is explicitly broken by 
gauge-invariant renormalizable terms — Yukawa cou- 
plings between fermion multiplets and the Higgs doublet 
and Majorana mass terms of the fermion singlets. The 
flavor symmetry-breaking Lagrangian is given by 

C = (Yu).. QjH - D? (Yd),. %i? t 



-E? (Y E \. LjH^ - N't (Y v )^ LjH 



--NfM u N5 + h.c, 



(6) 



where H = (1, 2)i/ 2 is the Higgs doublet, and gauge and 
Lorentz indices have been suppressed. The Yukawa cou- 
plings Yjj^d,e are n g x n g matrices, whereas the neutrino 
Yukawa coupling Y v is an n' g x n g matrix. The singlet 
neutrino Majorana mass matrix M is a symmetric n' g x n' g 
matrix. In the Standard Model without neutrino singlets, 
renormalizable terms proportional to Y v and M are ab- 
sent. 

Under the chiral flavor symmetry transformations ip — > 
U$ ip, where are unitary matrices in flavor space for 
the fermion fields ijj = Q, U c , D c , L, E c and N c , the 
Yukawa coupling matrices, the Majorana mass matrix 
and the angles transform as 

T v 7/ 



r MU N c, 

2 arg det — arg det Uu<= — argdetZ/foc, 

3 arg det Uq — arg det U l , 



Yu - 




Y D - 


■+ u D . 


Y E - 


-+ U E ' 


Y v - 


-> Un 


M - 


-> Un 


$3 - 


- 03 


$2 - 


- $2 


$1 - 


-> 01 



- arg det Ul — argdet£/_Ec. 



- arg det Ud" 
(7) 



Under CP, each matrix is transformed to its complex 
conjugate, and each •& angle changes sign, 



Yu,D,E,v 

M 

$1,2,3 



YjJ,D,E,v 1 

M* , 
-01,2,3 ■ 



(8) 



Under the chiral flavor symmetry transformation, the 
i9 angles are shifted by Eq. (7). The invariant angle 0qcd 
is defined by 



0qcd = $3 + arg det Yjj + arg det Y D 



(9) 



The analogous angles 0i, 2 can not be separately defined, 
but one can define an invariant ^-parameter in the elec- 
troweak sector 







EW 



0a + 20i 
+2 arg det Yg 



- arg det Y v + - arg det Y D 



(10) 
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After electroweak symmetry breaking, the QED i9-angle 
is 2 if qed = $ew- The factor of two arises because the 
generators for a non-abelian gauge theory are normalized 
to Tr T a T b = S ab /2. 

In the absence of electroweak symmetry breaking, 
there are n' g massive Majorana neutrino singlets with 
masses of O(M), the heavy Majorana neutrino mass 
scale, and all other fermions are strictly massless. It is 
natural that M be of order the GUT scale, the scale 
at which the GUT gauge symmetry breaks to the Stan- 
dard Model gauge group, under which the N c fields are 
uncharged. When the Higgs field gets a vacuum expec- 
tation value v/y/2, the Yukawa matrices generate Dirac 
mass matrices for the quarks and leptons, 

v 

m U.D.E.v = Yu,D,E,v — y= ■ (11) 

with the same flavor transformation properties as the 
Yukawa couplings. The Dirac and Majorana mass matri- 
ces of the (n g + n'g) left-handed neutrino fields combine 
to form a neutrino mass term 

~ST X {M H ) XJ Mj, l<l,J<n g + n' g (12) 

where the (n g + n' g ) x (n g + n'g) neutrino mass matrix 
Mj\f is equal to the symmetric matrix 

* ■ ( I ^ ) • < 13 > 

The (n g + n' g ) neutrino fields M% are (i/i,Nf). The 
(n g + n' g ) mass eigenstates of Eq. (13) give the Majorana 
mass-eigenstate neutrino fields, which are linear combi- 
nations of Ui and Nf. The heavy neutrinos with masses 
O(M) are predominantly N c with an 0(v/M) admixture 
of v, and the light neutrinos with masses 0(v 2 /M) are 
predominantly v with an 0(v/M) admixture of N c . 

A low-energy effective field theory can be obtained 
from the seesaw theory by integrating out the n' g heavy 
Majorana neutrino mass eigenstates. In this low-energy 
theory, the leading flavor symmetry-breaking Lagrangian 
is given by 

£ EFT = -Ui {Yu) tJ Q 3 H - Df (Yb) y QjH^ 
-Et {Y E ).. I.JI 

+ ^(L i H)(C 5 ) ij (L j H) + h.c, (14) 

where the coefficient of the dimension-five operator [9] is 
given by 

C 5 =YjM- 1 Y v (15) 

to lowest order in the 1/M expansion. When the elec- 
troweak gauge symmetry breaks, the dimension-five op- 
erator yields an effective n g x n g Majorana mass matrix 

m 5 = -C 5 v 2 /2 (16) 



for the (primarily) weak doublet neutrinos. Under the 
flavor symmetries and CP, the flavor matrices Yu,d,e 
and •& angles $1,2,3 of the low-energy effective theory 
transform under chiral flavor symmetry and CP as in 
Eq. (7) and Eq. (8), respectively, whereas C5 transforms 
as 

C5 — ► UjJ C5 Ul, 

c 5 - c; , (17) 

respectively. 

We will analyze the flavor structure of both the seesaw 
theory and its low-energy effective theory. The analysis 
depends only on the flavor transformation properties of 
the Yukawa coupling and Majorana mass matrices (i.e. 
the fermion mass matrices). Thus, it applies to any the- 
ory which has Dirac and Majorana mass matrices with 
the same transformation properties as given here, regard- 
less of whether the Dirac mass terms are proportional 
to Yukawa couplings in the theory, or are generated by 
some mechanism from more fundamental parameters of 
the theory. 

III. MASSES, MIXING ANGLES AND PHASES 

In this section, we define the mass and mixing param- 
eters of the high-energy seesaw theory and its low-energy 
effective theory. Most of the section is a review of well- 
known results, and serves to define the parameters and 
notation which are needed later. The mass matrices of 
the high and low energy theories in the weak eigenstatc 
basis are transformed to the mass cigenstate basis by fla- 
vor rotations to obtain the fermion masses and mixing 
matrices. The counting of mixing angles and phases for 
the case n' g = n g follows the analysis of Ref. [15]. The 
counting of physical parameters is given here for the cases 
n' g > n g and n' g < n g , for completeness. 

Any complex matrix M can be written in the form 
M = U A U' where U and U' are unitary matrices, and A 
is a diagonal matrix with real, non-negative entries. If M 
is also a symmetric matrix, then it can be written in the 
form M = M T = U T A U, where U is a unitary matrix. 

A. High-Energy Theory 

The flavor matrices of the high-energy seesaw theory 
are written in Eq. (6) in the weak cigenstate basis. These 
flavor matrices are related to the mass eigenstate basis 
by 



Y v = 




Ac/ Uu, 


Y D = 




Au Ud, 


Y E = 






Y v = 


Uatc 


A, U, y , 


M = 




T An U^. 
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where Ajjd.b, A„ and Ajv arc n g x n„ 



n' g x n g and 



diagonal matrices respectively, with real, non- 



negative entries; 
unitary matrices 



and Uatc 



and Uu.d.e.l* are n g 
and U 



N' 



are n g x n ff uni- 



tary matrices, which transform the mass eigenstate ba- 
sis to the weak eigenstate basis. Performing the chi- 
ral flavor transformation Eq. (7) with Uu<= T = \Ju ~ , 
U D , T = U D -\ U E J = U E '\ U q = Uu~\ U L = \J E ~\ 



and U 



A" 



U'. 



brings the flavor matrices to the form 



Yu = 


Ac/, 


Y D = 




Y E = 


A E , 


Y v = 


W' 1 A v V, 


M = 


A N , 


Ut/UrT 


-\ V = U 



(19) 



where Vckm = UyUrT 1 , V = \J„l) E and W = 
Ujv c_1 (U^yc) T are the three unitary matrices which de- 
scribe flavor mixing in the seesaw theory. Vckm is the 
Cabibbo-Kobayashi-Maskawa mixing matrix in the quark 
sector. As is well-known, this n g x n g matrix corresponds 
to the mismatch between the unitary field redefinitions 
on U and D in the quark doublets Q required to diagonal- 
ize Yu and Yd • V is the analogue of the CKM matrix in 
the lepton sector; it is the n g x n g matrix corresponding 
to the mismatch between the unitary field redefinitions 
on v and E in the lepton doublets L required to diago- 
nalize Y v and Ye- W is an n' g x n' g mixing matrix in the 
lepton sector corresponding to the mismatch between the 
unitary field redefinitions on N c required to diagonalizc 
M and Y„. 

To proceed further, it is necessary to consider the three 
cases n' g = n g , n' g < n g and n' g > n g individually. We 
first specialize to the case n' g = n g considered previously 
in Rcf. [15] and review the analysis given there. The anal- 
ysis is then generalized to the cases n' g ^ n g . The quark 
sector only depends on the number of quark generations 
n g , but the lepton sector analysis depends on whether 



"3i 



n g < n g 



or n^ > n g . 



The real diagonal matrices Au,d,e,v,n are invariant un- 
der the rephasings, 



A^ 
A„ 
A N 



-;<I> 



e "* A^ e 
e"**" A, e 1 
Vn A n n N , 



(<!»„ 



<P = U,D,E, 



(20) 



where &u,d,e,v are real diagonal matrices, and 77 at is a 
diagonal matrix with allowed eigenvalues ±1. Only ±1 
rephasings arc allowed for the Majorana fields N c . Under 
these rephasings, the mixing matrices Vckm, V and W 
transform as 



Vckm 



Vckm e 



i'l'i 



Matrices 


Masses 


Angles 




Phases 




Ac; 


n g 












A_D 


n g 












Vckm 





\n g (n g - 


1) \{n a 


-i)K 


-2) 


Total 


2n g 


2 n 9( n 9 ~ 


1) i(n 9 




-2) 



TABLE I: Parameters in the quark sector for n g generations. 
The Ay and Ad rows give the parameters if Yjj or Yn are 
considered separately, and the third row gives the additional 
parameters if both Yjj and Yd are considered together. There 
are (n a — l) 2 mixing parameters (angles plus phases), and a 
parameters. 



. 1) 
total of (n g 2 + 1) 



V 

w 



V e 



i>l>£ 



(21) 



Quark Sector: The parameter counting in the quark sec- 
tor is well-known, and is summarized here for complete- 
ness. The matrices Au and A^> each contain n g eigen- 
values, which correspond to the [/-quark and D-quark 
masses, respectively, and are CP even. The quark mix- 



ing matrix Vckm is an n g 



unitary matrix with n g 



parameters. It is conventional to divide these parame- 
ters into angles and phases — angles are even under CP, 
whereas phases are odd under CP. If the Vckm matrix 
is CP invariant, it is an n g x n g real orthogonal matrix 
with n g (n g — 1)/2 parameters. The unitary matrix Vckm 
has n g (n g — l)/2 angles and n g (n g + l)/2 phases, and can 
be parametrized by 



(22) 



where x is an overall phase, <E> = diag(0, 4>2, ■ • • , 4>n g ), 
and \& = diag(0, -02, • ' ' ? V'ng)- The phase redefinitions 
<J>;y and $£i of Vckm m Eq. (21) can be chosen to remove 
the 2n g — 1 phases x, 4>i, ipu * = 2, • • • , n g . 3 Thus, Vckm 
has n g (n g + l)/2 — (2n g — 1) = (n g — l)(n g — 2)/2 net 
phases. This counting of parameters is summarized in 
Table I. 

We choose a parameterization Vckm = V(6i,5i) in 
terms of a standard functional form V, where the n g (n g — 
l)/2 angles 0j G [0, 7r/2] and the (n g — l)(n g — 2)/2 phases 
Si E [0, 2tt). The CKM matrix for n g = 3 is given by [29] 



V(0 



12, Cl3> ^23 

ci 3 si 3 e 
1 





C23 S23 
-S23 c 23 



C13 



C12 S12 
-S12 C12 
1 



(23) 



There are n g phases each in $(7 and &d, but the transformation 
<E>[7 = &d <x 1 leaves Vckm invariant. 



6 



where Sj = sinf?i and Cj = cos^. It is now conventional 
to call the angles #23, #13, #12 rather than 61,92,63- The 
standard form Eq. (23) has det V = 1. 

Lepton Sector: The matrices An and Ae each have n g 
eigenvalues which are CP even. The lepton mixing ma- 
trices V and W are n g x n g unitary matrices, which can 
be parametrized by 

V = e l * e l * V{6i,6i) e**/ 2 , 
W = e tx ' e 1 *' V(6^) e 1 *'/ 2 . (24) 

We use the same standard functional form V as for the 
quark sector, but with different numerical values for the 
arguments 8i and Si. 4 The factor of two in VP and V?' will 
be explained below. 

The rephasing transformations $„, <S>e and t]n of 
Eq. (21) can be used to (i) eliminate x, x' an d ipi : (ii) re- 
strict ^ to the range [0, 2ir) rather than [0, 47r), and (iii) 
eliminate either $ or <&', &wi noi feot/i. It is convenient to 
use the same domain [0, 2tt) for all phases, which is why 
iff' was scaled by a factor of 2. 

First consider amplitudes which depend only on Y v 
and Ye, but not on M. In this case, the mixing matrix 
W is no longer observable and can be set to unity. The 
mixing matrix V has (2n g — 1) allowed phase redefini- 
tions: n from <& v , n from and minus one, because 
$ v = $e oc 1 does not change V. Thus, the parame- 
ter counting for the mixing matrix V is identical to that 
for Vqkm m the quark sector, with n g (n g — l)/2 angles, 
and (n g — l)(n g — 2)/2 phases. Similarly, for amplitudes 
depending only on M and Y v and not on Ye, the mix- 
ing matrix V is no longer observable and can be set to 
unity. The mixing matrix W has n g allowed phase re- 
definitions $„. Thus, there are n g (n g — l)/2 angles and 
n g (n g + l)/2 — n g = n g (n g — l)/2 phases. If the three 
matrices M, Y v and Ye are considered together, then 
the mixing matrices V and W together can have 2n g 
allowed phase redefinitions due to <&„ and As com- 
pared with the case of only V or only W, where there 
were 2n g — 1 + n g phase redefinitions possible, we have 
(n g — 1) fewer phase redefinitions, and hence (n g — 1) 
additional observable phases. These (n g — 1) additional 
phases occur because the same phase redefinition $„ was 
present for both V and W, and so cannot be chosen to 
remove phases from both V and W. Thus, there are an 
additional (n g — 1) phases if all three mass matrices are 
considered together. These phases can be included in ei- 
ther V or W . The standard form of the mixing matrices 
which uses the <fr u phases to eliminate the 4> phases from 
V is given by 

V = V(6i,Si), 

W = e~ l * V(C<50 e**'/ 2 , (25) 



4 The use of the same symbols 6i for the quark and lepton sectors 
should cause no confusion, since we do not need to deal with 
mixing in both sectors simultaneously. 



Matrices 


Masses 


Angles 


Phases 


Ajv 


n g 








A* 


Tin 








A B 


n g 








V : Y„,Y E 





\n a (n g - 1) 




W :M,Y V 





\n g {n g - 1) 


\n g (n g - 1) 


$ qt 1 








n g - 1 


Total 


3n 9 


n 3 {n g - 1) 


n 3 (n g - 1) 



TABLE II: Parameters in the lepton sector for n' g = n g gen- 
erations. The An, Av and A_b rows give the parameters if M 
or Y v or Ye are considered separately. The V and W rows 
give the additional parameters if both Y u and Ye, or both 
M and Y v are considered together, respectively. The last 
row gives the additional parameters to those in the previous 
rows when all three matrices M, Y v and Ye are considered 
together. There are 2n g (n g — 1) mixing parameters (angles 
and phases), and a total of n g (2n g + 1) parameters. 



whereas the standard form of the mixing matrices which 
uses the $„ phases to eliminate the $' phases from W is 
given by 

V = e^ V(6i,6i), 

W = V(0J,$ e 1 ^' 2 . (26) 

In Eq. (25), V has the canonical CKM form with n g (n g — 
l)/2 angles &i and (n g — l)(n g — 2)/2 phases Si, whereas in 
Eq. (26), W has the canonical PMNS form with n g (n g — 
l)/2 angles 6\ and n g (n g — l)/2 phases consisting of the 
(ri g — l)(n g — 2)/2 phases Si and the (n g — 1) phases 
ip'i. In either basis, there are (n g — 1) additional phases 
$ = $ — which cannot be removed, and are observable. 
This parameter counting for n' g = n g is summarized in 
Table II. 

$ Angles: Once the mixing matrices have been put in 
standard form, one can perform additional phase rota- 
tions which leave the mixing matrices invariant to elim- 
inate i9 angles. The only allowed transformation is an 
overall phase rotation with &u = <f>£> = <f>Q 1, i.e. baryon 
number. Under this phase transformation, 

^3 -► 03, 

$2 — > ^2 - 3n g <i>Q, 

"di -» 0i + \n g cf> Q . (27) 

The transformation leaves 03 and 02 + 20i unchanged, 
so there are two physical angles remaining: 0qcd, the 
strong interaction CP-angle in the basis where the quark 
mass matrices are real and diagonal, and $ew = 02+201, 
the electroweak CP-angle in the basis where the quark 
and charged lepton mass matrices are real and diagonal. 
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Matrices 


Masses 


Angles 


Phases 


A N 
A„ 

A B 

V : Y V ,Y E 
W :M,Y U 
$ gt 1 


n 9 
n 9 
n„ 










\n g {n g - 1) 


\{n g - n' g ){n g - n s - 1) 






\n g (n g - 1) — n'g + 1 

§<« - 1) 
<-l 

i(n 9 -n' g )(n g -n' g + l) 


Total 


n g + 2n' g 


TlfjTlg fig 


Tig Tig Tig 



TABLE III: Parameters in the lepton sector for n g fermion generations and n' g < n g neutrino singlets. The total number of 
parameters is equal to the sum of the first six rows minus the last row. The parameters in U n _„< are removed from V . 



Matrices 


Masses 




Angles 


Phases 


A N 


n' g 










A„ 


n g 










A E 


n g 










V : Y V ,Y E 







\n g {n g - 1) 


|(n a -l)(n 9 -2) 


W :M,Y U 







\n'g{n' g - 1) 


\n' g {n' g + 1) - n g 


i gt 1 










Tl ff - 1 


Un' g -ng 





IK 


- n g )(n' g - n g -l) 


|(n' 9 -n g )(n s - n g + 1) 


Total 


2n g + n' g 




n g Tl g Tl g 


Tig Tig — Tig 



TABLE IV: Parameters in the lepton sector for n g fermion generations and n' g > n g neutrino singlets. The total number of 
parameters is equal to the sum of the first six rows minus the last row. The parameters in U n i _ n are removed from W . 



2. n'g < n g 



For n' < rt g , the n' x n g diagonal matrix A^ can be 



"g ~~ '"gt u±i<= 11 g "3 

written as 



A,. 



A, 



where denotes the n' g x (n g — n' g ) zero matrix, and A^ 
is a diagonal n' g x n' g matrix with n' g real non-negative 
eigenvalues. This matrix is invariant under 



A„ 



K o 



e**" 



(29) 



where U ng -n' g denotes an arbitrary (n g — n' g ) x (n g — n' g ) 
unitary matrix. The rephasing transformations of the 
lepton mixing matrices are 



The additional unitary transformation matrix in 
Eq. (30) can be used to eliminate parameters in V. 
The parameter counting for n' g < n g is summarized 
in Table III. The number of CP-even parameters is 
(n g n' g + n g +n' g ) and the number of CP-odd parameters 



(28) is ( n g n 'g ~~ n g)i consistent with the results of Ref. [10]. 



For n' > n g , the n' x n g diagonal matrix A u can be 



9 

written as 



A„ = 



A„ 




(31) 



where denotes the (n' g —n g ) x n g zero matrix, and A v is a 
diagonal n g x n g matrix with n g real positive eigenvalues. 
This matrix is invariant under 



V 



w 



e" 1 *- 
EC 1 —/ 





"a/ 




V e 8 * E , 













"a/ 


Un' g -n g _ 








(30) 



instead of Eq. (21). 



(32) 

where U n ^- ng denotes an arbitrary (n' g — n g ) x (n' g — 
n g ) unitary matrix. The rephasing transformation of the 
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lcpton mixing matrices is 



V 

w 







^nj-n. 



W 7?JV. 



(33) 



instead of Eq. (21). 

The additional unitary transformation matrix in 
Eq. (33) can be used to eliminate parameters in W. 



The parameter counting for 
in Table 
(n g n' a + ? 



> n. 



is summarized 



is (n 



IV. The number of CP-even parameters is 
g + n' g ) and the number of CP-odd parameters 



e'-g 



n 



ah 



consistent with the results of Ref. [10]. 



Matrices 


Masses 


Angles 




Phases 




A E 


n g 












A 5 


n g 












fpMNS 





\n g {n a - 


1) 


2 n g( n g ~ 


1) 


Total 


2n g 


2~n g {n g — 


1) 


2 n a{ n a ~ 


1) 



TABLE V: Parameters in the lepton sector of the low-energy 
effective theory for n g generations. The Ae and As rows give 
the parameters if or ms are considered separately. The 
c^pmns row gives the mixing angles and phases of the PMNS 
mixing matrix. 



B. Low-Energy Effective Theory 

The flavor matrices in the low-energy effective the- 
ory are written in Eq. (14) in the weak eigenstate basis. 
These matrices are related to the mass eigenstate basis 
by 

Y v = U V c A v U v , 
Y D = U D c A D U D , 
Y e = U sc A E Us, 
C 5 = U'l A 5 U' v . (34) 

Performing chiral flavor transformations in the low- 
energy theory with Uu<= T = U[/c _1 , W D c T = Ud^ 1 , 
U E c = Ubc- 1 , U Q = U v -\ U L = Ue' 1 brings the 
flavor matrices to the form 

Yu = Ay, 

Yd = A D Vq KU , 

Y E = A E , 

C5 = (Ppmns) A5 U PMm , (35) 

where Vckm = Uc/Ud -1 and £7p mns = U'^U^ -1 are the 
two unitary matrices which describe flavor mixing in the 
low-energy effective theory. Vckm is the CKM mixing 
matrix in the quark sector. Ppmns is the PMNS mixing 
matrix in the lepton sector, which is the lepton mixing 
matrix which is physically measurable at low energies. 

The real diagonal matrices Au^e>,e,5 are invariant un- 
der the rephasings 

A 4 , -> e" 1 *"' A 4 , e***, ip = U, D, E, 

A 5 -> r\ v A 5 r) v , (36) 

which correspond to arbitrary phase redefinitions of the 
fermion mass eigenstate fields U c , D c , E c , U, D and 
E, and the discrete rephasings v — ► r\ v v, where r\ v is a 
diagonal matrix with allowed eigenvalues ±1 for the low- 
energy Majorana neutrino fields. Under these rephasings, 
the mixing matrices of the effective theory transform as 

Vckm -> e~^' Vckm e**«, 
Ppmns e~ l * E C/pmns ■ (37) 



The quark mixing matrix Vckm has the angles and 
phases given in Table I as before. The counting of pa- 
rameters in the lepton sector is summarized in Table V, 
and is well-known. C/pmns contains n g (n g — l)/2 angles 
61 . The number of phases of Ppmns is n g (n g + 1)/2 minus 
the n g phase redefinitions for a total of n g (n g — l)/2 
phases consisting of (n g — l)(n g — 2)/2 phases Si and 
(n g — 1) phases ipi. The canonical parametrization of 
C/pmns is 

f/pMNS = V(0i,6i) e 8 */ 2 , (38) 

* = diag(0,^ 2 , ■ • -An)- 

For n g = 3, the low-energy lepton mixing matrix is 
given by 

l/pmns = v(e?\e?\e?\5M) 
(10 \ 

x / a , (39) 

V e**™' 2 j 

where the superscript (U) denotes quantities in the 
PMNS matrix. 



IV. INVARIANT THEORY 

In the previous sections, we have discussed the param- 
eters (masses, angles and phases) for the low- and high- 
energy theories. We would like to analyze the theories 
using invariant quantities written directly in terms of the 
original parameters of the theory, the matrices Yu,d,e,v 
and M. The structure of the invariants is highly non- 
trivial, and depends in an interesting way on the number 
of generations. 

To study the invariants, it is useful to introduce several 
mathematical results from invariant theory [25-27]. The 
general problem is the following: one has a set of variables 
xi, . . . x n which transform (reducibly or irreducibly) un- 
der the action of a group G. The set of polynomials in 
{xj} with complex coefficients form a ring C[xi, . . . ,x n ]. 
The polynomial ring C[xi, . . . , x n \ is a free ring on the 
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generators x%, . . . ,x n , i.e. it is given by taking linear com- 
binations of all possible products of powers of the gener- 
ators with coefficients in C, and there are no non-trivial 
relations among the generators. 

The ring C[x±, . . . ,x n ] G C C[xi, . . . ,x n ) is the set of 
G-invariant polynomials, i.e. those polynomials which are 
unchanged by the action of G. This is clearly a ring, since 
sums and products of invariant polynomials are also in- 
variant polynomials. A highly non-trivial result, if G is 



a reductive group, is that C[xi, . 



is finite gener- 



ated. Let the generators be ii, . . . I r , each of which is a G- 
invariant polynomial in the original variables xi, . . . ,x n . 
Then, any G-invariant polynomial can be written as a 
polynomial P e C[/i, . . . , I r }. However, C[xi, . . . , x n ] G 
need not be a free ring in the generators I\ , . . . I r ; there 
can be non-trivial relations among them. 

In the following sections, we analyze the invariant ring 
for the quark and lepton sectors of the Standard Model 
effective theory and the seesaw model. It is useful to first 
look at some simple examples before discussing the case 
of interest. We start with a famous result on symmetric 
polynomials, and then discuss three examples involving 
continuous groups which arc closer in structure to the 
quark and lepton invariant problem. The first model is 
a theory which has a freely generated ring, with no rela- 
tions. The second theory has one non-trivial relation, and 
is similar in structure to the ring for quark invariants for 
three generations studied in Sec. VB and for lepton in- 
variants in the Standard Model for two generations stud- 
ied in Sec. VIA. The third example is only slightly more 
complicated, but leads to an intricate structure of in- 
variants, with many relations, and a complicated Hilbcrt 
series. This is similar to what we find for lepton invari- 
ants in the Standard Model for three generations, and in 
the seesaw model for two and three generations. 



A. Symmetric Polynomials 

The classic example from invariant theory is the study 
of symmetric polynomials. The permutation group S n 
acts on a polynomial f(x\, . . . , x n ) in C[xi, . . . , x n ] by 



p : f(xi, ...,x n )-> f(x p{1} , 



c p(n); 



(40) 



where (p(l), . . . ,p(n)) is a permutation of (1, . . . , n). A 
polynomial in C[xi, . . . ,x n ] Sn is invariant under the ac- 
tion of any permutation. A standard result [30] is that 
the invariant ring is generated by the elementary sym- 
metric polynomials 



h = Xx 



x 2 



5 A reductive group is defined by the property that every repre- 
sentation is completely reducible. A Lie group which is a direct 
product of simple compact Lie groups and U(l) factors is reduc- 
tive, as is any finite group. 



h = XiX 2 + X±X 3 + . . . + X n -\X n = ^2 Xi 
h = XiX 2 X 3 + ... 



i<3 



n — c lXn—\Xri ^ XiXjXf^ 1 

i<j<k 



X\X 2 ...X n 



(41) 



In other words, any symmetric polynomial f{x\, . . . , x n ) 
can be written as a polynomial in 1\ , . . . , /„ , 
f{xi,...,x n ) = g(h, . . . ,/„), e.g. 



2 I 2 , 2 

x 1 -\- x 2 -t . ■ ■ x n 



It - 27, 



(42) 



The important point is that g(Ii, ■ ■ ■ ,I n ) is a polyno- 
mial — otherwise the result would be trivial, for know- 
ing Ii, . . . , /„, one could solve Eq. (41) to find Xi, . . . , x n , 
and hence determine /. 



B. Model I 

Consider a theory with two couplings m\ and m 2 which 
transform under a G — U(l) x U(l) symmetry as 



m 1 



m 2 



e° vz m 2 . 



(43) 



We look at the ring C[mi, mj, m 2 , m 2 ] 
polynomials which are U(l) x U(l) invariant. It is clear 
that they can be written as linear combinations of mono- 
mials of the form 



(mimj) ri (m 2 m2) T 



(44) 



where r± and r 2 are integers. Thus, the ring of invariant 
polynomials is generated by the invariants I\ = m\m\ 
and I 2 = m 2 m2, and there are no relations between these 
generators. 

The Hilbert series H (q) is defined as 



H(q) = Y J C r q r 



(45) 



r=0 



where c r is the number of invariants of degree r, and 
Co = 1. In our example, c\ — 0; c 2 = 2 since m\m\ 
and m 2 m 2 are the two degree- two invariants; c 3 = 0; 
C4 = 3 since (miraj) 2 , (mim\)(m 2 m,2) and (m 2 m2) are 
the three degree-four invariants; and so on. It is easy to 
sec that the Hilbert scries is 

.6 , 8 , 



H(q) = l + 2^ + 3 9 4 + 4 9 b + 5g 

OO 

= ^(n + l) 9 2n 



n=0 



1 



(1-9 2 ) 2 



(46) 



Another derivation of the Hilbert scries is the follow- 
ing. The generators I\ = m\m\ and I 2 = m 2 m,2 are 
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both of degree two, and the invariants of higher order 
are given by multiplying together arbitrary powers of I\ 
and J2. The product 

(1 + h + It + ■ ■ •) (1 + h + I\ + ■ ■ ■) (47) 
gives each invariant once, which leads to the Hilbert series 



H(q) = {l + q' + q i + ...){l + q< + q l + 
1 



(l-q 2 ) 2 ' 



(48) 



in agreement with Eq. (46). 

In the general case of a scmisimple Lie group, it is 
known that H (q) has the rational form 



H(q) = 



N(q) 
D(q) 



(49) 



where the numerator N(q) and denominator D(q) are 
polynomials. Furthermore, the numerator is of degree 
(In and is of the form 



N(q) = 1 + cyq + . . . c dN ^q dN - X + q d > 



(50) 



where the coefficients are non-negative, c r > 0, and N(q) 
is palindromic, i.e. 



N(q) = q dN N(l/q). 
The denominator takes the form 



(51) 



(52) 



and is of degree do = d r . The number of denomina- 
tor factors p is equal to the number of parameters. The 
number of parameters is defined as the minimal codi- 
mension of an orbit, and agrees with the usual physics 
usage of the term. Model I has p = 2 parameters, be- 
cause we start with four objects TOi, m,2, to* and m*, 
(or equivalcntly, the real and imaginary parts of mi and 
TO2), and have two phase redefinitions Eq. (43), which 
eliminates two variables. In other words, one can always 
make a phase redefinition to make mi and m.2 real and 
non-negative, and these are the two independent param- 
eters. In our example, N(q) = 1, d\ = efe = 2 and the 
number of denominator factors is two. The number of 
denominator factors p is equal to the number of param- 
eters. 

There is a theorem due to Knop [31] which says that 



dim V > do — djv > 



V 



(53) 



where dim V is the dimension of the vector space on 
which the group transformations act; do and djv are the 
degrees of the denominator and numerator; and p is the 
number of parameters. In most cases, the upper bound 
is an equality, but not always. (We will see an example 
for the quark invariants involving only the [/-quark mass 



matrix.) In Model I, the vector space basis is mi, to*, TO2, 
TOj, so dim V = 4, p = 2, d^ = and do = d r = 4, 
and we see that Knop's theorem gives 4 > 4 — > 2, with 
an equality for the upper bound. 

One also can construct a multi-graded Hilbert series. 
Let c rir . 2r3ri be the number of invariants of order n in 
mi, order r-i in m*, order r% in 7712, and order in TO2 • 
Then 

92, 93,94) = 'Yl l Crir 2 r 3 r i q\ 1 ql 2 ql s ql i 
1 



(l-9i92)(l-9394)' 



(54) 



and the usual Hilbert series is H{q) = h(q,q,q,q). The 
multi-graded series gives more information about the 
structure of the invariants. However, it is important 
to remember that the results quoted above for H(q), 
Eqs. (49) (53), do not hold in general for the multi- 
graded case. 



C. Model II 

Consider a theory with couplings mi and 7712 with 
charges one and two, respectively, under a G — U(l) 
symmetry, 



TOi 



"mi, 



TO 2 



TO2 . 



(55) 



The ring of invariant polynomials C[mi, m\, TO2 , m^\ u ^ 
is generated by the four basic invariants I\ = toito^ , I2 = 
TO2TO2, I3 = m2m 1 2 and I4 = ra^m\. These generators, 
however, are not all independent, since 23/4 = I^Is, so 
that C[toi, to*, TO2, to*,]^*- 1 -* is not a free ring generated 
by h-4. 

It is straightforward to show that the multi-graded 
Hilbert scries is 



^(91,92,93,94) = 



1 



9i9l9394 



(1 - 9i92)(l - 9394)(1 - 9395)(1 



949?)' 
(56) 



where q\ , qi , qs and 94 count powers of mi , to* , TO2 and 
TO2, respectively. 

The denominator of the multi-graded Hilbert series is 
generated by the invariants ii_4, whereas the numerator 
compensates for the fact that I3I4 and l\l2 count as 
only one invariant at order q\q\qzq^ rather than two, 
because I3I4 = I\l2- The numerator of the multi-graded 
Hilbert series does not have the special properties of the 
numerator of the Hilbert scries H(q) discussed in the 
previous example. 

In this example, dimy = 4, dimG = 1, and there are 
three parameters since the phase transformation Eq. (55) 
eliminates one of the original four real variables in mi and 
TO2- The Hilbert series H(q) = h(q,q,q,q) is 



H(q) 



1 



9 



(l- g 2)2 (1 _ 5 3)' 



(57) 
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which has a palindromic numerator with djv = 3, and 
a denominator with do = 7, and p = 3 is equal to the 
number of denominator factors and to the number of pa- 
rameters. Knop's theorem gives 4 > 7 — 3 > 3, with an 
equality for the upper bound. 

Expanding Eq. (57) in a series in q gives the invariants 
of each degree. We see that there are two generators of 
degree two, I\ and h, and one generator of degree three, 
which can be chosen to be 7 3 + 7 4 , corresponding to the 
denominator factors (1 — q 2 ) 2 and (1 — q 3 ), respectively. 
Expanding out the denominator would give a coefficient 
of q 3 of +1. There are two invariants of degree three, 
I3 ± I4. The missing degree-three invariant 73 — 7 4 is 
counted by the +q 3 term in the numerator, so that the 
coefficient of q 3 in the expansion of H(q) is 2. When the 
denominator factors are expanded in a series, they can 
occur to any power, so one can have arbitrary powers of 
h, I2 and 13+14. However, the q 3 factor in the numerator 
occurs only once. This means that powers of -Z3 — 7 4 
higher than the first can all be eliminated in terms of 
polynomials P(h, h, I3 + I a) which have already been 
included. This statement follows from the identity 

{h-hf = [h+hf-Ahh 

= (h + hf - 4/?/ 2 . (58) 

There exists a similar identity for the Jarlskog invariant 
which will be derived in Sec. V. 

The generator 73 + 7 4 of the denominator is not ho- 
mogeneous in the multi-grading; 73 is of degree q^q 2 and 
7 4 is of degree q^q\, which is why Eq. (56) can not be 
written in a form similar to Eq. (57) with positive coe- 
ficients in the numerator and one less generator in the 
denominator. 



J 


™3 *2 

= m 2 m 3 , 


T 

h 




ho 


* 

= TO1TO2TO3, 


hi 


= m*m2?7J3, 


hi 


= mim^m^ 2 


ha 


* * 2 
= TO 1 m 3 TO 2 . 



There are 35 relations between products of invariants 
hlj given by: 7 4 7 5 = 7?7 2 , 7 4 7 7 = l\hi, 7 4 7 8 = I 2 I 2 , 
hh = I\2i hho = hh, hhi = hh2, lihz = hhho, 
hh = Iiho> hh = ^i3) ^5^9 = hliii hho = hh3, 
hhi = hh, h-hi = hhhi, hh — I 3 h, hh = Ifo? 
hh = hhh2, hhi = hhh, hh2 — hh 2 - hh3 = 
hl 2 , hh = hhhs, hh - I?!, hho - hhh, 
hhi = ^1^11, ^7^13 = hh 2 , hh = h\l 2 , hhi = hh-hz, 
hh2 = I 2 hhoi hho = hhh2, hh3 = h^hhi, 
ho hi = hhhi hoh2 = hhh, hoh3 = hh, hih2 = 
hhh, hih3 = hhh and 7 12 7 13 = 7 1 7|7 3 . The new 
feature here is that these relations are not independent — 
there are relations among the relations (known as syzy- 
gies in the mathematics literature), e.g. multiplying both 
sides of I4I7 — I 2 hi and 7s7g = 7^7io gives 

7 4 7 5 7 6 7 7 = 7^7 10 7 n , (61) 

which is also obtained by multiplying the relations 7 4 7s = 
I 2 I 2 and 7 6 7 7 = 7f7 3 , and using 7 10 7 n = hhh- The 
Hilbcrt series is 

1 + g 2 + 3g 3 + 4 g 4 + 4g 5 + 4g 6 + 3g 7 + g s + g 10 
~ (l-q 2 ) 2 (l-q 3 )(l-q*)(l-qS) 

(62) 



D. Model III 

Consider yet another model with three couplings mi, 
?7i2 and m 3 with charges 1, 2 and 3, respectively, under 
a U(l) symmetry, 

mi -> e^mi, m 2 -> e 2 ^m 2 , m 3 -> e 3i<#, m 3 . (59) 

The structure of the invariants is considerably more com- 
plicated than in the previous examples, even though the 
theory is only slightly more complicated. All the invari- 
ant polynomials are generated by thirteen invariant gen- 
erators 

7i = mim*, 
h = m2inl, 
7 3 = to 3 777,3, 

T 2 * 

I4 = 777 1 777 2 , 

T *2 

h = m\ml, 
h = m* 3 m 3 , 



Here dim\^ = 6, dimG = 1, and the number of param- 
eters is 5. From the Hilbert series, <1n = 10, do = 16, 
and p = 5. The number of parameters is equal to p, and 
Knop's theorem gives 6 > 16 — 10 > 5, with an equality 
for the upper bound. 

There are thirteen invariants in Eq. (60). However, 
there are only five denominator factors in Eq. (62), so 
only five basic invariants, two of degree two, and one each 
of degrees three, four and five, generate a free ring. The 
other invariants must satisfy non-trivial relations (those 
given below Eq. (60)), and this is reflected by the com- 
plicated numerator in Eq. (62), which implies that the 
invariant ring has a non-trivial structure, with many re- 
lations. The different terms in the numerator show that 
there are many invariants which can be eliminated when 
raised to higher powers, or multiplied by lower order in- 
variants, by relations analogous to Eq. (58). There is one 
invariant of degree two (the +q 2 term), three in degree 
three (the +3g 3 term), etc. This model shows that even 
a relatively simple theory given by Eq. (59) can lead to 
a set of invariants with an interesting syzygy structure. 
Furthermore, the number of invariants and relations of 
each degree is encoded in the Hilbert series. 
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V. QUARK INVARIANTS 

We can now address the first problem of interest - 
flavor invariants in the quark sector. We are interested 
in polynomials in mu, mu^ , and where 



mu 
m D 



Ujjc mu Uq, 
U D o T m D Uq, 



(63) 



under the chiral flavor transformations. 6 To cancel Uu<= 
and l/flc, one must consider the combinations 



Xjj = mu^mu, 
X D = niD^mD, 

which both transform as adjoints 

Xu,d — > Uq Xu,D Uq. 



(64) 



(65) 



Thus, the invariants are traces of products of Xjj and 
Xjj. The structure of the invariants depends non- 
trivially on the number of generations, so we consider 
the cases n g = 2 and n g = 3 separately. 



First, consider invariants involving only Xjj . The basic 
invariants are 



(Xu),(Xu 2 ) 



(Xu 3 ) 



(66) 



where (*) denotes a matrix trace. This series of traces 
terminates after n g terms for an n g x n g matrix, by the 
Cayley-Hamilton theorem which states that every matrix 
satisfies its characteristic equation. For an arbitrary 2x2 
matrix A, the Cayley-Hamilton theorem gives 



A 2 = (A) A 



(a 2 ) (Ay 



1. 



(67) 



Taking the trace of both sides gives the trivial result 
(A 2 ) = (A 2 ). Multiplying by A and taking the trace 
implies 



(A 3 ) = 1(A) (A 2 ) 



\(Af 



(68) 



so that (A n ), n > 3 can be written in terms of (A) 
and (A 2 ). Thus, there are two independent invariants, 
-^2,0 = (Xu) an d Ia.q = (Xi/ 2 ), which can be constructed 
from Xjj alone. Both of these invariants are CP even. 
The two invariants contain the same information as the 
eigenvalues of Xjj, i.e. the two [/-type quark masses. For 
invariants constructed only from mu, the number of pa- 
rameters is p = 2, the two eigenvalues of Xu- The vector 



One could equally well work with the Yukawa matrices, which 
differ by factor v / v2. 



space has dimU = 8, because mu and mu^ are both 
2x2 matrices, and ^2,0 and -£4,0 are of degree two and 
four, respectively, in mu, so the Hilbert series is 



H(q) 



(l-q 2 )(l-q±)- 



(69) 



Here djj = 0, dp = 6 are the degrees of the numerator 
and denominator, respectively, and the number of de- 
nominator factors is p = 2, which is equal to the number 
of parameters. Knop's theorem gives 8 > 6—0 > 2, which 
holds, but this time the upper bound is not an equality. 

Similarly, there are two independent CP-even invari- 
ants Jo. 2 = (Xd) and Jo, 4 = (Xd 2 ) which involve only 
Xjj. These two invariants contain the same information 
as the eigenvalues of Xjj, namely the two -D-type quark 
masses. 

Invariants containing both Xu and Xjj can be written 
as traces of the form 



(Xu ri X D s 'Xu r2 X D s 



(70) 



for integers r» and s,. The Cayley-Hamilton theorem for 
a 2 x 2 matrix, Eq. (67), implies that all powers Ti and 
Si greater than one in Eq. (70) can be reduced, so we are 
left with traces of the form 



(XuXjj ■ ■ ■ XuXjj) — ((XuXjj) v ) 



(71) 



Again, invariants with r > 1 can be rewritten in terms of 
lower order invariants, so there is only one independent 
invariant, 7 2 ,2 = (XuXjj), which is CP even. 

In summary, the basic quark invariants for n g = 2 
quark generations, which generate all the invariants, are: 

h,o = (Xu) = (mu^mu) , 

I0.2 = (X D ) = (m D ^m D ) , 

h,o = (Xu 2 ) = ((mr/Vy) 2 ) , 

^2.2 = (XuX D ) = (mu ] mum D ] m D ) , 

io,4 = (Xd 2 ) = ((m D ^m D ) 2 ) . (72) 

Writing the invariants in terms of the usual quark masses 
and the Cabibbo angle gives 

h,a = m 2 u + ml, 
Io,2 = ml + m 2 ., 

T 4,4 

h,2 = mlml + m\m\ + ( m * s - m^m 2 . - m* u ) cos 2 9, 

(73) 



T 4,4 

J ,4 = m d + m s 



Knowing the five invariants allows one to determine the 
four masses and 9, because mi > 0, and 9 lies in the first 
quadrant. 

Using u and d to count powers of mu and mD gives 
the multi-graded Hilbert series 



h(u, d) 



1 



(1 - u 2 )(l - u 4 )(l - d 2 ){\ - rf 4 )(l - u 2 d 2 ) ' 



(74) 
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The Hilbert series H(q) 
H(q) = 



h(q,q) is 
1 



(l- 9 a)2(l_g4)3- 



(75) 



In this example, p = 5 (four masses and one mixing an- 
gle, see Table I), dimV = 16, since there arc four 2x2 
matrices, djv = 0, and do = 16. The number of denom- 
inator factors is the number of parameters, and Knop's 
theorem gives 16 > 16 — > 5, with the upper bound an 
equality. 

The denominator factors in Eq. (75) show that there 
are two generators of degree two, and three of degree 
four, which agrees with Eq. (72). 

If one started with Xjj and Xjj as the basic objects, 
then dim V = 8. In this case, the Hilbert series is given by 
replacing q 2 — > q in Eq. (75), since we now count powers 
of Xjj, Xjj rather than mjj,mjj, so = 0, djj = 8 and 
Knop's inequality becomes 8 > 8 — > 5. 



For arbitrary 3x3 matrices A, B and C, one has the 
identity 







(A) 2 (B) (C)-(BC) (A) 2 — 2 (AB) (A) (C) 



-2 (AC) (A) (B) + 2 (ABC) (A) + 2 (ACB) (A) 
- (A 2 ) (B) (C) + 2 (AB) (AC) + (A 2 ) (BC) 
+2 (C) (A 2 B) + 2 (B) (A 2 C) - 2 (A 2 BC) 

-2 (A 2 CB) - 2 (ABAC) (78) 
which can be derived by substituting A — > A + B + C into 
Eq. (77), and picking out the order A 2 BC terms. This 
identity eliminates (ABAC), i.e. traces where the same 
matrix is repeated, so that in invariants Eq. (70), Xjj, 
Xjj 2 , Xd and Xjj 2 can each occur at most once. For 
example, (Xjj . . . Xjj . . .) can be replaced by (Xjj . . .), 
and (Xfj . . . Xfj . . .) can be replaced by (Xfj . . .), which 
can then be eliminated using Eq. (76). 

Writing out all of the possibilities gives the basic quark 
invariants for n g = 3 quark generations. There are 11 
CP-even invariants, ten of which arc 



For an arbitrary 3x3 matrix A, the Cayley-Hamilton 
theorem states that 



A'' 



A 2 (A) 



1 



(A) 2 (A 2 ) 



l r 
6 



(A) 3 - 3 (A 2 ) (A) + 2 (A 3 ) 



1. 



(76) 



Taking the trace of both sides gives the trivial result 
(A 3 ) = (A 3 ). Multiplying by A and taking the trace 
gives 

(A 4 ) = \(At-(A) 2 {A 2 ) + A -{A 3 )(A) + \{A 2 ) 2 , 

(77) 

so that (A n ), n > 4 can be rewritten in terms of (A), 
(A 2 ), and (A 3 ). 

Thus, the invariants involving Xjj alone are I2.0 — 
(Xjj), /40 = (Xjj 2 ) and Igfi = (Xjj 3 ), and invariants 
involving Xjj alone arc /0.2 = (Xd), Ioa = (Xd 2 ) and 
Io,e = (Xd 3 ), all of which are CP even. 

Invariants containing both Xjj and Xjj are of the form 
Eq. (70), but now with n = 1,2 and s, = 1,2, so that 
one has traces of products of Xjj, X^, Xo, X\\. This re- 
striction still leads to an infinite number of invariants. 
However, many of these invariants are not independent. 



-^2,0 


= (Xu), 


-^0,2 


= (Xd), 


-^4,0 


= (Xu 2 ), 


-^2,2 


= (XjjX d ) , 


-^0,4 


= (Xd 2 ), 


-^6,0 


= (Xu 3 ), 


h,2 


= (Xu 2 X D ) , 


h,4 


= (XuX D 2 ) , 


Io,6 


= (Xd 3 ), 


h,4 


= (Xjj 2 X d 2 ) , 



(79) 



and one CP-odd invariant 

I e q = (Xu 2 Xd 2 XjjXd) — (Xe> 2 Xjj 2 XdXjj) 



The eleventh CP-even invariant is 

^6^5' = (Xu 2 Xd 2 XjjXd) + (Xd 2 Xjj 2 XdXjj) 



(80) 



(81) 



All the invariants in the quark sector can be written as 
polynomials in these basic invariants. 

The multi-graded and one-variable Hilbert series are 



h(u, d) 



H(q) - h(q,q) 



1 + u e d 6 

(1 - U 2 ){1 - M 4 )(l - W 6 )(l - d 2 )(l - d 4 )(l - d G )(l - U 2 d 2 )(l - U 4 fi 2 )(l - U 2 (i 4 )(l - U 4 d 4 ) ' 

1 + q 12 



(l-q 2 ) 2 (l-q±) 3 (l-qSni-q*y 



(82) 



14 



respectively. This case has p = 10 parameters, consisting 
of 6 masses, three angles and one phase, which agrees 
with the number of denominator factors. The original 
variable space has dimV = 36, from the two 3x3 mass 
matrices and their complex conjugates. The degrees of 
the numerator and denominator are d^ = 12 and d]j = 
48, respectively, and Knop's inequality is 36 > 48 — 12 > 
10, which is satisfied, with the upper bound being an 
equality. If one started with Xjj and Xn as the basic 
objects, then dimU = 18, and the Hilbert series is given 
by replacing q 2 — > q in Eq. (82), so rfjv = 6, do = 24, 
and Knop's inequality becomes 18 > 24 — 6 > 10. 

The denominator of Eq. (82) shows that there are two 
invariants of degree two, three of degree four, four of 
degree six, and one of degree eight, which can occur 
multiplied in arbitrary combinations, with no relations 
among them. This is a total of 10 invariants, which are 
the ones listed in Eq. (79). One can see that their de- 
grees match the denominator factors in Eq. (82). What 
about the remaining two invariants? The numerator fac- 
tor of Eq. (82) shows that there is one additional in- 
variant of degree twelve other than those given by prod- 
ucts of denominator factors. This is the CP-odd invari- 
ant Eq. (80). The Hilbert series implies that the other 
degree-twelve invariant, Eq. (81), cannot be an indepen- 
dent invariant. Indeed, it can be written as a polynomial 
in the other CP-even invariants, 

= ^2,0^0,2 _ ^2,0^4,0-^0,2 _ 3/2,2-^2,0^0,2 

+3/4,2/2,0-^0,2 — -^0,4/2,0-^0,2 + 3/2,4/2,0-^0,2 
— 3/4,4/2,o/o,2 

+-^0,4/6,0-^0,2 + 3/2,4/4,2 + 3/2,2/4,4 

+/o, 6/2, 0/4,0 — /o, 6/6,0, (83) 

and so can be eliminated. 

The Hilbert series numerator only has an entry q 12 , but 
there is no q 2i term. This means that /g 6 is an indepen- 
dent invariant, but the square and all higher powers of 
/g g'' arc not. The square of the CP-odd invariant /g § is 
CP-even, and can be written as a polynomial (with 241 
terms out of a possible 305 terms) in the CP-even invari- 
ants of Eq. (79). The most general polynomial invariant 
in the quark sector can be written as 

P1+1UP2 (84) 

where Pi and P2 are polynomials in the CP-even invari- 
ants Eq. (79). 

This example illustrates how the structure of the in- 
variants is encoded in the Hilbert series. For many pur- 
poses, the details of the relations, such as Eq. (83), or the 

formula for ^/g 6 ^ are not important; all one needs to 

know is that /g 6 ^ occurs linearly, and /g g can be elimi- 
nated. 

The quark sector parameters are determined by the 
ten CP-even parameters / 2 ,o, /4,o, /6,o, /o,2, /o,4, /o,6, 



/2.2 ; /2.4, /4,2, Z4.4, and the single CP-odd parameter 
/g 6 '. From the CP-even invariants, one can determine 
the [/-type quark masses m u ^ c _t and -D-type quark masses 
md,s,b, which are real and non-negative, and four combi- 
nations of the CKM parameters, cos 6*12, cos 6*13, cos #23 
and cos S, all of which are CP even. Since the CKM an- 
gles 0i2, 013, 023 lie in the first quadrant, these angles 
arc determined uniquely by their cosines. However, cos S 
does not determine the phase 5 uniquely, because it can- 
not distinguish between S and —S. Under CP, 5 <-> —8. 
Thus, one Z2 piece of information, the sign of 5, is miss- 
ing. This sign is provided by the invariant /g 6 ^ . The 
only information needed is the sign of /g 6 ^ , which is why 

^/g 6 ^ can be written in terms of the other CP-even in- 
variants. This discussion corresponds to the well-known 
result that the unitarity triangle can be obtained by mea- 
suring the lengths of its sides, which are CP-conserving, 
rather than the angles, which are CP-violating. Know- 
ing the sides determines the triangle up to a two-fold 
reflection ambiguity, which is fixed by the sign of Iq~q , 
or, equivalently, the sign of the Jarlskog invariant, so that 
the only additional information contained in the Jarlskog 
invariant is the sign. The relations between the invari- 
ants are similar to those obtained by studying rephasing 
invariants [15]. 

The invariant /g 6 also can be written as 

1$ = l([Xu,X D } 3 ) , (85) 

and is proportional to the Jarlskog invariant J [12], 

4" } = 2iJ(m 2 -m 2 u ){m 2 -m 2 c ){m 2 -m 2 u ) 

x(m 2 s -m 2 d )(m 2 b -m 2 s )(m 2 b -m 2 d ), (86) 

where 

,/ = Im (Vckm)h {Vckm)* 12 (^ckm) 2 2 (^ckm)2i ■ (87) 

/g 6 ^ vanishes if two U -type quarks or two -D-typc quarks 
are degenerate. It is well-known that quark CP viola- 
tion vanishes for degenerate [/-type or D-type quarks. 
iitf is odd under the exchange of two ?7-type or two D- 
type masses, e.g under m u <-> m c , whereas the invariants 
in Eq. (79) are even under exchange, so /g 6 ^ cannot be 

written in terms of the other invariants, ^/g 6 ^J is even 

under exchange, and can be written in terms of the other 
invariants. 

It is, of course, well-known that CP conservation in 
the quark sector requires J = 0, or equivalently, /g 6 = 
0. What is new is the structure of the ring of all in- 
variant polynomials, and the relation between the CP- 
conscrving and CP- violating invariants. 



VI. LEPTON INVARIANTS FOR TWO 
GENERATIONS 

The structure of the lepton invariants, like the quark 
invariants, depends on the number of generations, so we 
first consider the case of n g = 2 generations in this sec- 
tion. The case of n g — 3 generations is considered in 
Section VII. We will outline the derivation of the results, 
but not give all the details. 

A. The Standard Model Effective Theory 

We now study the lepton invariants in the Standard 
Model low-energy effective theory with a neutrino Ma- 
jorana mass term. The structure of the lepton invari- 
ants is considerably more complicated than the quark 
invariants. The lepton sector of the low-energy theory 
contains the flavor symmetry breaking matrices Y E and 
C5, so we are interested in polynomials in ttie, tre* , 771,5 
and TO5* = 777,5 since 7775 is a symmetric matrix. These 
matrices transform as 

tob — > U E <= T m E U L , 

777 E t — > Ue^ m E ^ U L * , 

m 5 — > U L T m 5 U L , 
m 5 * Uj} 7775* U L * , (88) 

under chiral flavor transformations. To cancel Ue", one 
must consider the combinations 

X E = m E ^m E , 
X E = X E T = m E T m E *, (89) 

which transform as 

X E -> U{ X E U l , 
X e t -» U L T X E T U L * . (90) 

It also is convenient to define 

X 5 = 7775*7775, (91) 

which transforms as 

X 5 -> Ui) X 5 U L , (92) 

as well as ^7775* (X E n ) T 7775^ , which transforms as 

(7775* (X E n f m B ) -> (7775* (X E n ) T m 5 ) U L . 

(93) 
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The invariants involving only Xe are ^2,0 = {Xe) and 
74 = (X E 2 ), whereas the invariants involving only 7775 

and 7775* are Jo, 2 = (-^s) and 7o,4 = (^5 2 )- 

The invariants involving Xe, 7775 and 7775* are of the 
form 

(7775* (X E ri f m 5 X E sl ...7775* (X E r "f m 5 X E s -) 

(94) 

for integers r% and s%. The Cayley- Hamilton theorem im- 
plies that all powers fj and s% greater than one in Eq. (94) 
can be rewritten in terms of lower order invariants. Thus, 
one needs to consider traces of matrix products contain- 
ing the matrices Xe, X$, and (7775* X E T 7775) at most 
once. 

In summary, the generators of the invariants are: 

72,o = {X E ) = (m E ^m E ) , 
7o,2 = {X 5 ) = (7775*7775) , 

7 4:0 = (X E 2 ) = {(m E ^m E ) ) , 

7 2 ,2 = (m 5 * X E T 777 5 ) = (TO5 X E 777 5 *) 

= (m E T m E * 7775 7775*) = (777 B t 777_b 7775* 7775) , 
7o,4 = (X 5 2 ) = ((7775*7775)') , 

74,2 = ("75* X E T 777 5 X E ) 

= (7775* m E T m E * 777 5 ms'TriE) , (95) 

4^4 = ( m 5* X E T 7775 X E 7775* 7775) 

- (7775* X E T 7775 7775* 7775 X E ) 

= (7775* m E T m E * 7775 ms'TriE m 5 * m 5 ) 

- (7775* m E T m E * 7775 7775* 7775 m E 'm E ) , 

where 7j ± is CP odd, and the rest are CP even. The 

square of the CP-odd invariant, (i^^ , is not inde- 
pendent; it can be expressed in terms of polynomials in 
the other CP-even invariants. In addition, the CP-even 
invariant 7^, obtained by the substitution ► + in 

l[ 4 , is not independent, and thus is not included in the 
above list. 

There are six parameters: four masses, one angle and 
one phase, see Table V. The four masses, one mixing 
angle, and one phase, can be determined from 72,0, li.o, 
7o,2; 7o,4, 72.2 and 72,4 up to a sign ambiguity in the phase, 
just as for the case of three generations of quarks already 
discussed. The sign of the phase is fixed by the sign of 
7 M 

The multi-graded Hilbert series is 



h{y,z) 



1 + 7/z 4 

(1 - y 2 )(l - ?/)(! - z 2 )(l - z 4 )(l - y 2 z 2 )(l - y*z*) ' 



(96) 
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where y counts powers of and z counts powers of 777,5. 
The single variable Hilbert series is 



H(q) = h(q,q) = 



1 



(I_g2)2(l_g4)3(l_g6) 



(97) 



The q 8 term in the numerator shows that there is one 
degree-eight invariant 4 which occurs, but that the 
square of this invariant is not independent and can be 
eliminated. 

The number of denominator factors p = 6 is equal to 
the number of parameters, and djv = 8, dp = 22. The 
number of variables is dim V = 14, since we have one 2x2 
mass matrix, one 2x2 symmetric mass matrix, and their 
complex conjugates. Knop's inequality 14 > 22 — 8 > 6 
is satisfied, with an equality for the upper bound. The 
six parameters correspond to 2 charged lepton masses, 
2 Majorana neutrino masses, one mixing angle and one 
phase. 

The denominator of Eq. (97) shows that there are two 
generators of degree two, three of degree four, and one of 
degree six, which agrees with the CP-even invariants in 
Eq. (96). The numerator shows that there is an invariant 
of degree eight, whose square can be eliminated, which is 
lj± J . The structure of the invariants for n g = 2 is similar 
to that for quarks for n g = 3. 

Weak-basis invariants for two generations in the 
low-energy effective theory were studied previously by 
Branco, Lavoura and Rcbelo [32]. They defined an in- 
variant Q, related to l[ J by 

2iImTrQ = 1$ , (98) 

and showed that Q = is a necessary and sufficient con- 
dition for CP conservation. This is consistent with our 
results, since the only CP-odd generating invariant is 

1 IA ■ 



B. The Seesaw Model 



In this section, we analyze the lepton invariants in the 



seesaw theory for n fl 



2 generations of fermions. 



There are three matrices in the lepton sector, m„, mg 
and M, and their complex conjugates m v \ rtiE^ and 
M' = M*J From Eq. (7), we see that only 777 £• trans- 
forms under Ue", so it must always occur in the combi- 
nation 



X E = mE^rriE, 
which transforms as 

X E -» Ul X E U L 



(99) 



(100) 



7 It is worth emphasizing that in our notation m„ refers to the 
Dirac mass matrix m v = Y v vj \/2, not the Majorana mass matrix 
7715 of the effective theory. 



under the chiral flavor symmetry transformations. The 
mass matrices m v , ml, M and M* transform as 



ml 
M 
M* 

It is useful to define 

X„ 
Z, 



U N c 777„ U L , 

Ui) ml U N c*. 
U N c T MU N o, 

Uj^c^ M* lijqc 



(101) 



m' v m v , 

t 



v 777 ;y 777 ^ , 

Z v = Z v = m v m v 



(102) 



which transform as 



as well as 



Z N 
Zx 



which transform as 

Xn 
Zn 
Zx 



U{ X v U L , 

lAjqc Zj^ IAn° 



M*M, 
MM*, 
m,, X R mJ 



U Nc Xn Un<=, 
Un<? Zn IAn c 
Un<? Zx IAn°* 



(103) 



(104) 



(105) 



Note that Z N L = Z N * = X N . 

The invariants involve three mass matrices, m e, m v 
and M. One first can consider the simpler problem of 
studying invariants which only depend on two out of the 
three matrices. The first case, invariants involving only 
777^ and m u , consists of invariants formed from traces 
of Xe and X v only, with no insertions of M or M* . 
These invariants are the same as the invariants in the 
quark sector with the substitutions X\j X v and Xe — > 
Xe- The second case, invariants involving only m v and 
M, are invariants which do not contain Xe- These have 
the same structure as invariants constructed in the low- 
energy theory, with the replacements 7775 — > M, mE — > 
777^, i.e. X E — » Zj. 

The most general invariant involving all three matrices 
has the structure 



(M*A 1 MA 7 2 



.M*A 2n _ 1 MA^ n ) 



(106) 



where Aj = 1 or A, = m v T'{XE,X v )m 1 / , where V is 
a polynomial in Xe and X v . This result can be ob- 
tained by representing the chiral transformations of the 
matrices graphically, as shown in Fig. 1. Products of 
matrices such as Eq. (106) also occurred when studying 
rephasing invariants [15]. For rcphasing invariants, one 
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,t V- - 



M* 



FIG. 1: Graphical representation of the chiral transformation 
properties of the lepton mass matrices Xe, rn v , ml, M and 
M* . A solid line represents Un" , and a dashed line Ul- The 
invariants are obtained by forming graphs with no external 
lines. 



can factor long products into smaller ones, each involving 
at most four mixing matrices, using reconnection identi- 
ties. This factorization is no longer possible for the case 
of mass-matrix invariants, which leads to an interesting 
and highly non-trivial structure for the invariants. 

The basic invariants can be constructed using Eq. (106) 
and eliminating higher powers of matrices by the Cayley- 
Hamilton identity Eq. (67). The generators are: 



£2,0,0 
A), 2,0 
-£o,0,2 

-£4,0,0 
-£2,2,0 
-£0,4,0 
-£0,2,2 
-£0,0,4 
-£2,2,2 
-£0,4,2 
-£2,4,2 

t(-) 



0,4,4 

£4,4,2 
2,4,4 



(X E ) = (mE^mE) , 

(X v ) = (mjm v ) , 

(X N ) = (M*M) , 

(X E 2 ) = (rns t mBTns t ms) , 

(X V X E ) = (mJm u mE*m E ) , 

(X v 2 ) = {mj m u mj rriv) , 

(Z U Z N ) = (m.mjMM*) , 

(X N 2 ) = (M*MM*M) , 

(ZxZn) = (m^mE TUE^n-v MM*) , 

(M*Z V MZ V T ) = (M*m v mjMm v *m v T ) , 

(M*Z U MZ X T ) 

(M* m v m v ^ Mm v *mE T i r nE*'i r n u T ) , 
(M*Z U Z X M) - (M*Z X Z V M) 
(M* m v mj m^rriE^ TriEm,) M) 

— (M* m^mE^ mEm^ m u m^ M) , 
(Z N Z V MZ V T M*) - (M*Z v Z N MZj) 
(MM*m u mjMm l/ *m l/ T M*) 

- (M*m u mJ MM*Mm v *m v T ) , 
(M*Z X MZ X T ) 



(AT 



m v mE rnEfn 



THE ms m u ) 



2,6,2 



/ ( - ) 

4,4,4 

4,6,2 



(M M* rnvrriE^ ttieITIv ^Mm v *m v T M*) 

— (AI*m v mE • mETn^ MM' Mm v *m v T ) , 
(M*Z v Z x MZj) - (M*Z X Z V MZ V T ) 
(M* m v mj m v mE^ vriErnJ Mrn v * m v T ) 

— (M* 'm u mE ^mETfi^ ! m u m u ^ Mm v * 'm v T ) , 
(M*Z N Z x MZl) - (M*ZlZ N MZ x ) , 

(M*Z V Z X MZ'£) - (M*Z X Z V MZ T X ) . 



(107) 



There arc several invariants which can be immediately 
eliminated because they are polynomials in lower order 
invariants and which have not been listed above. These 

r(+) r(+) r(+) r(+) r(+) 
z 2,4,2' J 0,4,4' J 2,4,4i J 2,6,2' ^4,4,4 

j|^ 2 , which are related in an obvious way to the invari- 
ants in Eq. (107) with superscripts (— ). The degree-eight 



invariants £2^2 anc ^ -£q 4 



r(+) 



are eliminated by the identities 



— £o,0,2£o,2,o£2,0,0 — £o,0,2£o,4,0-£2,0,0 — 2/o,2,0-£2,2,2 



2/o,2,2-£2,2,0 + 2/ 2 ^ 2 , 

-'0,4,2 — i O,0,4-'o,2,0 



' 2/q, 0,2^ 



r2 r2 
J 0, 0,2^0, 2,0 

— 2£q,2,2 + 2£q^4j 



(108) 



and the degree-ten invariants I%a a anc ^ -£262 are elimi- 
nated by the identities 



— £o,0,2£o,2,o£2,2,0 ~ 2/0,0, 2£2, 4, 2 _ £o,0,4£o,2,0-£2,2,0 

— 2/o,2,2£2,2,2 + 2/ 2 ^ 4 , 
= £q,2 q£o,2,2£2,0,0 ^ 2/o,2,o£2,4,2 ~ £o,2,2£o,4,0-£2,0,0 



-2/o, 4,2-£2,2,0 + 2/ 2 ^ 2 - 



(109) 



(Z N Z X MZ V T M*) - (M*Z X Z N MZ U 1 



The degree-twelve invariants /j 4 4 and I4 g 2 are a ^ so 
polynomials in lower order invariants, but we do not in- 
clude the explicit identities here. 

In Eq. (107), there are three CP-even invariants of 
degree two, five of degree four, two of degree six, one of 
degree eight, and one of degree ten, for a grand total of 12 
basic CP-even invariants. In addition, there are two CP- 
odd invariants of degree eight, two of degree ten and two 
of degree twelve, for a total of 6 basic CP-odd invariants. 
All of the invariants can be written as polynomials in 
these 18 basic invariants. 

The multi-graded Hilbert scries is 



u \ N 

1 + 2x 2 y A z 2 + y 4 z 4 + x 2 y i z i + x 2 y 6 z 2 + x^y^z 1 + x' ± y v z A — x"y"z" — x~y"z* — x'y"z" ~ x'y"z' 



N 



„4„,4 4 



„4„,6„2 



„2„,6 „6 



„2„,8 _4 



„4„,6 fi 



„4„,8 _4 
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-x 6 y s z 4 - 2x 4 y s z 6 - x 6 y 12 z s , 
D = (1 - x 2 ) (1 - x 4 ) (1 - y 2 ) (1 - y 4 ) (1 - z 2 ) (1 - z 4 ) (1 - *V) (1 - l/V) (1 - z W) (l - y 4 z 2 ) (l - x 4 y 4 z 2 ) , 

(110) 



where cc, y, z count powers of the, m v and M, respec- 
tively. The Hilbcrt scries H(q) = h(q,q,q) is 

1 + g 6 + 3g 8 + 2g 10 + 3g 12 + g 14 + g 20 
{q) (l-g 2 ) 3 (l-g 4 ) 5 (l-g 6 )(l-g 10 ) ' 

(111) 

which has a palindromic numerator. The number of de- 
nominator factors p = 10 is equal to the number of pa- 
rameters, and d/v = 20 and <4> = 42. The number of 
variables is dimT^ = 22, because we have two 2x2 ma- 
trices with 4 independent entries, one 2x2 symmetric 
matrix with 3 independent entries, and their complex 
conjugates. Knop's inequality is 22 > 42 — 20 > 10, and 
the upper bound is an equality. The 10 parameters in 
the lepton sector of the seesaw model for n g = n' = 2 
generations correspond to 2 charged lepton masses, 4 Ma- 
jorana neutrino masses of the two light and the two heavy 
neutrinos, 2 angles and 2 phases. 

One can see from the Hilbert series that the struc- 
ture of invariants is far more complicated than in the 
quark case. The denominator factors (1 — g 2 ) 3 (l — g 4 ) 5 
of Eq. (Ill) corresponds to the generators 4,o,o, 4, 2,0, 
4,o,2, 4,o,o, 4,2,o, 4,4,o, 4,2,2, 4,o,4- At degree six, in 
addition to products of lower order invariants, there are 
two new invariants, 4.2,2 and 4.4.2- These two invari- 
ants correspond to the (1 — g 6 ) factor in the denomina- 
tor, and the +g 6 term in the numerator. Since there is 
only one power of (1 — g 6 ) factor in the denominator, we 
know that there will be non-trivial relations involving the 
degree-six invariants. At degree eight, there are 3 new in- 
variants from the +3g 8 term in the numerator in addition 
to products of lower degree invariants which make up the 
denominator. These are the three degree-eight invariants 
in Eq. (107). There are three new invariants of degree 
twelve (from the +3g 12 ), but only two degree-twelve in- 
variants in Eq. (107). The third degree-twelve invariant 
is the square of the degree-six invariant corresponding 
to the +g 6 term in the numerator, so the square of this 
CP-even invariant cannot be removed. We have noted 
earlier that there must be non-trivial relations involving 
the degree-six invariants. These relations first occur at 
degree 14, 

= 4,0,24,2,0-^2,6,2 + 4 2 ,2, 0^0,4,4-^2, 0,0 ~ ^0,2,0-^2,4,4 

— 4,2,o4, 2,24, 4,2 — 4, 2, o4,4,44, 2,0 — 24,2,2^,6^2 
— 4,4,04,4,44,0,0 + 24,4,04^,4,4 + 24,4,2^,4^2 

= 4 2 ,0,24,2, 0-^2,4^2 — 4 2 ,0,24,e!2 + 4,0,24,2,0^,4^4 

— 4,0,24,2,24^2 — 4, 0,24,2, o4,4,4 ~ 4,0,44,2,0^,4^2 
+ 2/o,0,442, 6^2 — 2/0,2,242,4^ + 24, 2, 24^4^4, (112) 



and are non-linear relations involving the two degree-six 
invariants. One can proceed to higher degrees — there 
are six relations of degree 16, etc., and verify the num- 
ber of independent invariants at each degree agrees with 
Eq. (111). The details of the relations are not important. 
The main purpose of giving Eq. (112) is to show that 
there can be non-linear relations among the generating 
invariants. To completely unravel all of the non-linear 
relations requires going beyond degree 20, the highest 
power of g in the numerator of Eq. (111). 

VII. LEPTON INVARIANTS FOR THREE 
GENERATIONS 

In this section, we consider the lepton invariants in 
the low-energy and high-energy theories for three gen- 
erations of fermions. The number of invariants is far 
greater than for two generations, and there are many re- 
lations between them. For the low-energy theory, we give 
the Hilbcrt series, and the invariants which correspond to 
the denominator factors. For three generations, even the 
Hilbert series proved too difficult to compute. For this 
case, we make some general remarks, and discuss some 
invariants considered previously by Branco et al. [32, 33], 
and by Davidson and Kitano [24]. 

A. The Standard Model Effective Theory 

The invariants involving only Xe are /2.0 == (Xe), 
4o = (Xe 2 ) and 4,o = (Ag 3 ), whereas the invariants 
involving only m.5 and 7715* are 4,2 = (A5), /o,4 = (A5 2 ) 
and 4, 6 = (X 5 3 ). 

The invariants involving Xe, and 7715* are of the 
form 

(m 5 * (X E ri f m 5 X E Sl -..m 5 * (X E r "f m 5 X E Sn ) 

(H3) 

for integers ri and Si. The Cayley-Hamilton theorem 
implies that all powers r t and Si greater than two in 
Eq. (113) can be rewritten in terms of lower order invari- 
ants. Thus, one needs to consider traces of matrix prod- 
ucts containing the matrices Xe, X 5 , (rn 5 * X E T rn^), 

and ^?7i5* (X E 2 ) T m^j at most twice. Identity Eq. (78) 
cannot be used to eliminate traces with multiple powers 
of m 5 , because (m^Am^B) gets converted to traces of the 
form (m 2 AB^ which are no longer invariant. There are 
many basic invariants, which involve a single trace, up 
to degree m^m 1 ^ , and we do not list them all here. The 
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ones up to degree twelve, which are sufficient for the de- 
nominator of the Hilbert series (and hence to determine 
the parameters) arc: 

^2,0 = (X E ) = {m E ^m E ) , 
Io,2 = (X 5 ) = (m 5 *m 5 ) , 

hfl = (X E 2 ) = {{mE^mE) 2 ) , 

h,2 = (X E X 5 ) = (m E ^m E m 5 *m 5 ) , 

I ,4 = (X 5 2 ) = ((m 5 *m 5 ) 2 ) , 

^6,0 = (X E 3 ) = {{m E ] m E f) , 

J 4,2 = (Xe 2 X 5 ) = ((TO £ t TO £ ) 2 m 5 *m 5 ) , 

h,2 = (m 5 * X E T m 5 X E ) 

= (m 5 * m E T m E * m 5 m E ] m E ) , 

h,i = (X E X 5 2 ) = (niE^niE (m 5 *m 5 ) 2 ) , 

Io,6 = (*5 3 ) = ((m 5 *m 5 ) 3 ) , 

h,2 = (m 5 * X E T m 5 X E 2 ) 

= (m 5 * m E T m E * m 5 (m E ^m E ) 2 ) , 

= (m 5 * X E T m 5 m 5 * m 5 X E ) 

I 



± (m 5 * m 5 m 5 * X E T m 5 X E ) 
= (7775* m E T mE* ^5 7715* 7775 77773^777/3} 
± (7775* 7775 m.5* "niE T mE* m§ m-E^mE) , 
ha = (m 5 * (X E T ) 2 m 5 X E 2 ) 

= (m 5 * (■m E T m E *) 2 m 5 (m B t m B ) 2 ) , 
4ji = (™5* X E T m 5 m 5 * m 5 X E 2 ) 

± (m 5 * 7715 m 5 * X E T m 5 X E 2 ) 
= (7775* mE T m E * m 5 7775* 7775 (mE^mE) 2 ) 
± (7775* 7775 777.5* m E^ m E* m-5 (77773^77173) ) , 
4a = (™5* (X E T ) 2 m 5 7775* 7775 X E 2 ) 

± (7775 * 777 5 777.5* (X B T ) 2 777 5 X B 2 ) 
= (t77 5 * {m E T m E *) 2 7775 »775* "75 ( 777 B t 777 B ) 2 ) 

±(7775*77757775* (m E T mE*) 2 m 5 (m^m^) 2 ) . 

(114) 

The multi-graded Hilbert series is 



= ^, 

iV = -y 24 z 18 - 2y 20 z 14 - 2t/ 20 z 12 - y 20 z w - 2t/ 8 z 14 - 3y 18 z 12 - y 18 z 10 - 3y 16 z 14 - 3y 16 z 12 - 3t/ 6 z 10 - y 16 z 8 
~y 16 z 6 - y 14 z 14 - y 14 z 12 - y 14 z 10 - 2y 14 z 8 - y 14 z 6 - y 12 z 14 + y 12 z 4 + y w z 12 + 2y 10 z 10 + y w z 8 + y 10 z 6 
+y w z 4 + y 8 z 12 + y 8 z w + 3y 8 z 8 + 3y 8 z 6 + 3y s z 4 + y 6 z 8 + 3y 6 z 6 + 2y 6 z 4 + y 4 z 8 + 2y 4 z 6 + 2y 4 z 4 + 1, 

D = (1 - y 2 ) (1 - y 4 ) (1 - 7/) (1 - z 2 ) (l - z 4 ) (l - z 6 ) (l - y 2 z 2 ) (l - y 4 z 2 ) 2 (l - y 2 z 4 ) (l - y e z 2 ) 

x (l - y 4 z 4 ) (l - y 8 z 2 ) , (115) 

where y counts powers of 77773 and z counts powers of 7775. The single- variable series H(q) = h(q, q) is 

_ 1 + g 6 + 2q 8 + 4g 10 + 8g 12 + 7q 14 + 9g 16 + 10g 18 + 9g 20 + 7g 22 + 8q 24 + 4g 26 + 2 g 28 + g 30 + q 36 

(l-g 2 ) 2 (l-g 4 ) 3 (l-g 6 ) 4 (l-g 8 ) 2 (l-9 10 ) ' 



The number of denominator factors p = 12 is equal to 
the number of parameters, and cIn = 36 and (i/j = 66. 
The number of variables is dim V = 30, because we have 
one 3x3 matrix with 9 independent entries, one 3x3 
symmetric matrix with 6 independent entries, and their 
complex conjugates. Knop's inequality is 30 > 66 — 36 > 
12, and the upper bound is an equality. Note that the 
numerator is palindromic. The 12 parameters consist of 3 
charged lepton masses, 3 Majorana light neutrino masses, 
3 angles and 3 phases. 

The Hilbert series Eq. (116) has a complicated numer- 
ator, which shows that the structure of the invariant ring 
is highly non-trivial. From the denominator of Eq. (116), 
we see that there are two generators of degree two, three 



of degree four, four of degree six, two of degree eight, and 
one of degree 10, which can be multiplied freely, with no 
relations. These account for most of the invariants in 
Eq. (114), but there remains one CP-even invariant each 
of degrees 6, 10, 12, and one CP-odd invariant each of 
degrees 8, 10, 12. These contribute q 6 + q 8 + 2q 10 + 2q 12 
to the numerator in Eq. (116). The coefficient of q 8 in 
the numerator of Eq. (116) is 2. Where does the other 
degree-eight invariant not in Eq. (114) come from? The 
degree-six invariant that corresponds to the numerator 
factor q e can be multiplied by cither of the two degree 
invariants, ^2,0 or Jo, 2, to give two additional degree-8 in- 
variants. One of these can be written as a polynomial in 
lower order invariants; the other survives. One can con- 
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tinue this analysis to arbitrarily high order — the entire 
invariant structure is encoded in a very compact way in 
the Hilbert series Eq. (116). An explicit example of the 
construction just discussed is given in Sec. VI B for the 
high-energy theory with n g = 2, which provides a simpler 
example of an invariant ring with non-trivial relations. 

For three generations, Branco, Lavoura and Rcbelo [32] 
defined four invariants: 

2ih = i$ 

2il 2 = i^XEm^m^m\m^rnt > X] s 'm<^ — c.c. 

21/3 = l^X ETn^m^m^m^m\X]^m^rn% l m^ — c.c. 

2H4 = det \m5XErri5 + 7715X^7715] — c.c. (117) 

of degrees (4,4), (4,6), (4,8) and (6,6), and showed that 
the vanishing of these invariants implies CP conserva- 
tion. The CP-violating invariants of Eq. (114) corre- 
spond to the denominator factors of the Hilbert series. 
There are additional CP-violating invariants not listed 
which correspond to terms in the numerator. 



B. The Seesaw Model 

The invariants involve three mass matrices, m,E, m v 
and M. One first can consider the simpler problem of 
studying invariants which only depend on two out of the 
three matrices. The first case, invariants involving only 
niE and m„, consists of invariants formed from traces 
of Xe and X v only, with no insertions of M or M*. 
These invariants are in direct analogy to the invariants 
of the quark sector with the substitutions Xjj — > X v 
and Xo — > Xe- The second case, invariants involving 
only m„ and M, are invariants which do not contain Xe- 
These have the same structure as invariants constructed 
in the low-energy theory, with the replacements 7775 — > 
M, the — > 777y , i.e. Xe — * . 

The most general invariant involving all three matrices 
has the structure 

(M*A X MA% . . . M*A 2n . 1 MA T 2n ) , (118) 

where Ai = 1 or Ai — ■m v V{XE-,X v )m u ^ . where V is 
a polynomial in Xe and X v . The generating invariants 
are given by using Eq. (118). In this case, there are a 
very large number of generating invariants. They include 
all those discussed earlier in the seesaw theory for two 
generations, as well as many other. 

For n g — n' g — 3 generations, there are 21 parameters 
which consist of 9 masses, 6 angles and 6 phases. The 9 
masses are the 3 charged lepton masses, 3 light Majorana 
neutrino masses and 3 heavy Majorana neutrino masses. 
There are 3 angles in the mixing matrix V and 3 angles 
in the mixing matrix W. There is one <5-type phase in V 
and in W, two Majorana phases W in W, and 2 phases $ 
which are not removeable when V and W are considered 
together. 



We have been unable to construct the multi-graded 
and one-variable Hilbert series in this case. However, it 
is clear that the structure of the invariant relations is ex- 
tremely complicated. There are a number of constraints 
on the form of the one-variable Hilbert series. The de- 
nominator must be a product of p = 21 factors. The 
numerator must be palindromic, and d/v and cLd must 
satisfy the Knop inequality 48 > do — e?jv > 21 since 
dim V = 48. The number of variables dim 1/ = 48 results 
because there are two 3x3 matrices the and m„ with 9 
independent entries each, one 3x3 symmetric matrix M 
with 6 independent entries, and the complex conjugates 
of the three matrices. 

Ref. [33] defined six invariants in the seesaw theory, 

2ih = (Y U Y}M*MM*{Y„Y}) T M) -c.c. 

2ih = (Y u Y 1 }M*MM*MM*(Y u Yj) T M) - c.c. 

2%h = (Y u YjM*MM*MM*(Y u Y2) T MM*M) -c.c. 

(119) 

which involve CP- violating phases which are relevant for 
leptogenesis, as well as 

2ih = (Y„X E YjM*MM*(Y v X E Yj) T M)- c.c. 

2ih = (Y U X E Y}M*MM*MM*{Y U X E Y^) T M) -c.c. 

2ih = (Y U X E Y}M*MM*MM*{Y U X E Y^) T MM*M) -c.c. 

(120) 

which involve the other phases. 
Ref. [24] defines an invariant 

2ih = (KWtyj^rv^r 1 ) (121) 

for leptogenesis, where k is 7775 with factors of the Higgs 
vacuum expectation value removed. This is not a poly- 
nomial in the basic variables of the seesaw model. It can 
be related to the invariants considered here using the for- 
mulas given below. 

Invariants in the seesaw model can be related to those 
of the low-energy effective theory. The basic relation is 
Eq. (15), which relates the neutrino mass matrices in 
the seesaw model to the Majorana mass matrix 7775 in 
the low-energy effective theory. Clearly, the relations be- 
tween the invariants cannot be polynomial, since inverse 
powers of M are involved, but one can write the low- 
energy invariants in terms of a rational function of the 
high-energy invariants. The basic identities are: 

det A A' 1 = (A)- A 

dot A = I(A) 2 -i(yl 2 ) (122) 

for 2x2 matrices, and 

det A A^ 1 = A 2 - A(A) -^{A 2 } + ^(A) 2 

dot A = l(A 3 )-±(A 2 )(A) + ±(A) 3 (123) 
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for 3x3 matrices, which can be combined with 

C 5 = YjM~ x Y v = Y?{M*M)- 1 M*Y V (124) 

to obtain the desired relations using A — M*M, and 
substituting for C5 (i.e. 777,5) m the expressions for the 
low-energy invariants. The expressions are valid as long 
as det M*M ^ 0, i.e. as long as the singlet neutrinos are 
heavy and the transition to a low-energy effective theory 
is valid. 



VIII. CONCLUSIONS 

We have used the mathematics of invariant theory to 
classify the independent invariants of the Standard Model 
effective theory and its high-energy seesaw model and 
to study the non-trivial structure of relations (syzygies) 
among the invariant generators. The complete classi- 
fication of invariants and the Hilbert series have been 
obtained for the Standard Model effective theory with 
a dimension-five Majorana neutrino mass operator. A 



complete solution also has been obtained for the renor- 
malizable seesaw model with n g = n' g = 2 fermion gen- 
erations. The lepton sector of the seesaw model involves 
three different mass matrices, the charged lepton mass 
matrix, the Dirac Mass matrix of the weakly-interacting 
doublet neutrinos and the Majorana mass matrix of the 
gauge-singlet neutrinos. The invariant structure is very 
complicated. In the case of n g = n' g = 3 generations 
of fermions, we have been unable to find the Hilbert se- 
ries for the invariant generators, and thus the structure 
of the syzygy relations for three generations remains an 
open problem. 
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